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I ■ ABSTRACT 

O \ We discuss the prospects for the detection of gamma-ray bursts (GRBs) by the 

Cherenkov Telescope Array (CTA), the next generation, ground-based facihty of imag- 
^ ■ ing atmospheric Cherenkov telescopes (lACTs) operating above a few tens of GeV. By 

' virtue of its fast slewing capabilities, the lower energy threshold compared to current 

lACTs, and the much larger effective area compared to satellite instruments, CTA 
CsJ ■ can measure the spectra and variability of GRBs with excellent photon statistics at 

\ multi-GeV energies, which would revolutionize our understanding of the physics of 

Op ■ GRBs, test their validity as the origin of ultra-high-energy cosmic rays, and provide 

\l ' powerful probes of the extragalactic background light as well as Lorentz-invariance 

violation. Employing a model of the GRB population whose properties are broadly 
consistent with observations by the Gamma-ray Burst Monitor (GBM) and Large 
\ Area Telescope (LAT) onboard Fermi, we simulate follow-up observations of GRBs 

with the Large Size Telescopes (LSTs), the component of CTA with the fastest slew 
speed and the best sensitivity at energies below a few hundred GeV. For our fiducial 
assumptions, we foresee that the LSTs can detect ~ 0.1 GRBs per year during the 
prompt phase and ~ 0.5 per year in the afterglow phase, considering only one array 
site and both GBM and the Space-based multi-band astronomical Variable Object 
Monitor [SVOM) as the alert instruments. The detection rates can be enhanced by 
^ ■ a factor of about 5 and 6 for the prompt emission and the afterglow, respectively, 

assuming two array sites with the same sensitivity and that the GBM localization 
error can be reduced to less than 1°. The expected distribution of redshift and photon 
counts are presented, showing that despite the modest event rate, hundreds or more 
multi-GeV photons can be anticipated from a single burst once they are detected. We 
also study how the detection rate depends on the intrinsic GRB properties and the 
delay time between the burst trigger and the follow-up observation. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are the most violent explo- 
sive phenomena in the Universe. The prompt emission of 
GRBs typically has luminosity of L-y ~ 10^^-10^^ erg s"^ 
at ~ 0.1-1 MeV, and the following long-lasting afterglow 
emission is observed in the radio to X-ray bands (see re- 
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views, iMeszarod [200^ : IZhand |2007^ . High-energy gamma- 
ray emission above ~ 0.1 GeV was first seen by the 
Energetic Gamma-Ray Experiment Telescope (EGRET) 
detector o n-board the Comp t on Gamma Ray Ob serva- 
tory (e.g., iHurlev et af] 1 1994 : iGonzalez et~all 120031 '). The 
Large Area Telescope (LAT) on-board the Fermi satellite 
has recently detected hi gh-energy gamma r ays f rom a much 
larger sample of GRBs "jAbdo et all I2009al lbll3. |2010| . I2OIII : 
lAckermann et all l20ld . 120111 ). The LAT GRBs exhibit the 
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following features (se e recen t revi e ws, iGranot et all I2OI0I : 
iBhat fc Guirie3 I2OIII : iPe'er I I2OIII : llnoue et al.l I2012I '). fa) 

Most of the LAT GRBs do not show significant suppres- 
sion at the high-energy end of their spectra, though ob- 
served LAT limits on the GeV fluence for the GRBs de- 
tected by the Gamma-ray Burst Monitor (GBM) onboard 
the Fermi satellite ma y suggest a steeping o r cutoff in the 
high-energy spectrum ( |Beniamini et aLll201 j ): (b) Some of 
the LAT GRBs h ave an anomalous extra component in the 
> 0.1 GeV range (|Abdo et al.ll20093 : lAckermann et al bold . 
[20l3), while others show a high-energy spect rum consistent 
with the Band function (|Band et al.l 1 199a ): (c) In many 
cases, the emission onset in the > 0.1 GeV energy range 
is delayed relative to that in the < 1 MeV energy range, 
(d) The LAT GRBs often show the long-lived high-energy 
emission lasting longer than the duration of the sub-MeV 
component; (e) Not only long GR Bs but also short GRB s 
seem to have the above features l|Ackermann et aLll2010l ). 
Understanding these features is likely to give us important 
clues to GRB mechanisms and related astrophysics. 

The Cherenkov Telescope Array (CTA), under plan 
as t he next generation internatio nal gamma-ray observa- 
tory ||CTA Consortiumll2010l . l201ll ). will provide a great step 
forward in studying these issues. CTA will be constructed 
at two sites, one each in the northern and southern hemi- 
spheres, and comprises three types of telescopes: the Large 
Size Telescopes (LSTs) with 23 m diameter and 4.6° field of 
view (FOV); the Medium Size Telescopes (MSTs) of 12 m 
and 8° FOV; and the Small Size Telescopes (SSTs) of 7 m 
and 10° FOV. With CTA, the sensitivity will improve by a 
factor of 5-10 in the 0.1-10 TeV range compared to existing 
Cherenkov telescopes. It will cover about 5 decades in en- 
ergy, including energies below a few tens of GeV and above 
100 TeV, and the angular and energy resolution will be ap- 
preciably increased. In the following, we summarize open 
issues that can be unraveled by CTA. We refer the reader to 
llnoue et al.l l|2012l ') for an extensive overview on the science 
prospects for GRB observations with CTA. 

(i) The bulk Lorentz factor V is one of the key quantities 
in understanding the properties of relativistic jets making 
GRBs. It is thought to be limited to order of ~ 1000 in 
the classical fireball model with baryons, whereas higher 
Lorentz f actors may be ach ieved by Poynting-dominated 
jets (e.g., ISpruit et al.l I2OOII ) or radiatio n-dominated jets 
with dissipation via e.g. , jet-confin ement l|loka et al.ll201ll ) 
or baryon-entrainment l|lokall2010l ). The bulk Lorentz fac- 
tor can be constrain ed from observations o f the high-energy 
end of the spectra (|Lithwick fc Sarill200ll '). and absence of 
high-energy cutoffs for the LAT GR Bs have indeed given us 
lower limits on V (e.g., Abdo et al. |[2009a,b,d). In the case 
of GRB 090926A (|Ackermann et al.ir2011i '). a sharp soften- 
ing of the spectrum is observed at ~ 1.4 GeV, and inter- 
preting this as the pair-creation cutoff leads to F ~ 720. 
However, such estimates c an be affected by the finite ex- 
tent of the emission region |Barind 20061). time-dependence 
of th e photon field (^Grano t et al.l l2008l : Asano fc IMesz arod 
2011 ) and multiple emis sion regions (|Li 2010; Aoi ct ~al] 



2010l : lHascott et al. Il2011^ ■ so that the resulting high-energy 



(ii) The mechanism of the prompt emission is one of the 
most critical issues in the theory of GRBs. In the classical 
scenario, ~ MeV emission is explained by the optically-thin 
synchro tron radiation from elec trons accelerated at internal 
shocks iRees fc Meszaro3ll994l ). However, this classical sce- 
nario h as several problems in expla ining observations (see 
reviews, iMeszarod^OO^ : IZhang||2007l ') , and many alternative 
scenarios s uch as the photosphcric dissipat ion scenario (e.g. 
Thompsonlll994l: iM eszar os fc Rocs 20qol: iRees fc Meszaroi 



20051 : lloka et al.l l2007l : iBeloborodovl I2OI0I : iMurase et al 



2012 ) and magnetic dissipation scenarios (e.g.. iLvutikov 
2006i : IZhang fc Yaiil l201lf : iMcKinnev fc Uzdenskvl I2OIID 



have been suggested. For the mechanisms of the > 0.1 GeV 
emission, both leptonic and hadronic emissio n have been 
have been propo s ed. Synchrotron em ission (|Wang et al.l 
I2OO9I : llokal bold : iDaigne et all bOld). synchrotron self- 



Compton (SSC) emission (jCorsi et al. |20 10a,b), exter- 
nal in verse-Compton (EIC) emission I Toma et a l| l200d . 
l2011bl ). an d proton synchrotron or proton-induced cascade 
emissions (lAsano et al1l2009bl . I2OI0I : iRazzaque et al.ll20ld : 
iMurase et al.ll2012l ) are currently viable explanations. How- 
ever, once the prompt emission is detected by CTA, it 
may be possible to discriminate among the emission mech- 
anisms by analyzing the very-high-energy (VHE) spectrum 
and variability with high photon statistics. 

(iii) It is believed that the long-lived high-energy compo- 
nent is related to the afterglow emission. In the simplest sce- 
nario, the high-energy emission is explained by synchrotron 
emission from the non-radiative external forward shock with 
extreme parameters (jKumar fc Barniol Duran 2011) or ra- 
diative external forward shock (jGhisellini et al.ll20ld ). The 
Klei n-Nishina effect m ay play a role in the initial rapid de- 
cay (|Wang et al.ll20l"ol ). but the more natura l explanation 
is the gradual turn-o f f of the prompt emiss ion (|Liu fc Wand 
, 201ll : iHe et akilioill : iMaxham et"alll201ll) . One of the im- 
portant tests for synchrotron external shock scenarios is to 
see the maximum sync hrotron cutoff that d ecreases as the 
Lorentz factor declines (jPiran fc Nakail20ld ) , and CTA may 
eventually see the SSC emission component. 

Observations by CTA can also provide a clue to 
the origin of X-ray shallow-decay emission. The shallow- 
decay behavior is observed i n m o st X-ray afterglows 
of Swift GRBs (iNousek et all l200d : lO'Brien et all I2OO6I : 
IZhang et al. I I2OO6'), though it was not predict e d by 
the standard a fterglow model (|Meszaros fc ReesI Il997i : 
ISari et all 1199^ ). Its origin has been debated for many 
years, and possible models include energy injection (e.g., 
Dai fcJLu 1998; Rccs & Mcszaros 1998; Zhang ot al. 200(|), 
the long-lasting internal activitv (|Ghisellini et al.l 120071 : 
iKumar et al. 20081: Murase et al] 201 ll), the long-lasting re - 



shock Icenet et al.ll2007l : lUhm fc Beloborodovll2007l ). 
time-dependent microphysics (*Io ka et al] 20061 ) . effect of 
cosmi c ray esca pe (Derrner 20q3)7 
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spectra can be more complicated. For more detailed studies, 
observations with much better statistics are required. CTA 
will be ideal for these purposes. 



sion (|Panaitescuil2008al ). prior e mission jY amazaki 200!^, 
and multi-compo nent jets (e.g., lEichler fc Granot, i2006l : 
iToma et al.l l200d ) . Th ese models shou l d predict different 
high-energy emissions fran et aL I l2008l : iPanaitescul l2008bl : 
[Murase et al. 2010, 2011j), so that the simultaneous observa- 
tions by Swift and the Space-based multi-band astro nomical 
Variable Object Monitor {SVOM) (|Paul et al.ll201ll ) as well 
as CTA would be useful for discriminating among the mod- 
els. On the other hand, high-energy emission has been ob- 
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served concurrently with X-ra y flares, which are often seen 
in the early afterglow phase l|Abdo et al.ll201lh . Although 
they are considered to originate from some internal dissiga- 
tion processes rather than external shocks (.loka et al.ii2005l ) , 
their mechanism and the origin of the high-energy emission 
is still unclear. CTA may provide a breakthrough for reveal- 
ing these problems. 

(iv) VHE signals from GRBs can be useful for other 
purposes. GRBs may be the main sources of the ob- 
served ultrahigh-e nergy cosmic rays (UHECRs) (jWaxmanl 
1 19951 : IVietrilll995l ). If this hypothesis is true, VHE gamma 
rays provide one of the crucial probes, and detections 
of char acteristic signals, including hadronic cascade ra- 
diation (iBottcher Dermeijll998l: (Pe'er fc WaxmanI l2005l: 



Dermer fc AtovanI l2006l : lAsano et al.l l2009al; [ Murase et al.l 
2012'), ion synchrotron radiation ( Bottcher fc Dermerlll99^ 



Totani 1998a, idi), synchrotron pair-echo em ission induced 
by ultra-high-energy photons (iMurasdbOll ). and emission 
via nucle ar photodisintegration an d Bethe-Heitler pair pro- 
duction l| Murase fc Beacomllioiol ). are relevant to test this 
GRB-UHECR hypothesis. Also, VHE gamma rays from 
GRBs are useful as a probe of the Universe. Since they 
are attenuated by the extragalactic background light (EBL), 
one can constrain the EBL by measuring the attenua- 
tion as in the cas e of blazars (e.g., lAharonian et al. I I2OO6I : 
[Albert et al.ll2007l ) . VHE gamma rays also lead to secondary 
emission as a result of interactions with the EBL, and de- 
tections or non-detections of this pair-echo emission en- 
able us to probe the uncertain interg alactic magnetic fields 
in voids as well as the EBL (e.g. , Razzague et al.l l2004l : 



iTakahashi et al.ll2008l : lMurase et aLll2009l ). GR Bs may occur 
even at the very distant un i verse, z > 5 (e.g., Kawai et al.l 
I2OO6I : ISalvaterra et al.ll2009l : iTanvir et all 120091) . and VHE 
gamma rays might be gen erated by GRB s originating from 
population III stars (Tom a et al.l 201 lal ). If we detect the 
attenuated VHE gamma rays from such distant GRBs, it 
would be possible to constrain the EBL in the high-redshift 
universe (|Gilmore et al.|[2009l : llnoue et al. 201G|) and /or cos - 
mic magnetic fields at high redshifts ( Takahashi et al.ll201lh . 
In addition, one can obtain insight into fundamental physics. 
For example, using the temporal (and spectral) information, 
one can test the Lorentz-invarianc e violation predicted in 
some quantum gravity theories fe.g..|Amelino-Camelia et al.l 
1 19981 : lAbdo et al.ll2009dl : IShao fc Mall2010l ). 



As described above, CTA should be a powerful tool to 
understand GRB physics and test other various ideas and 
hypotheses related to GRBs. VHE gamma-ray astronomy 
has now been firmly established by state-of-the-art lACTs 
such as H.E.S.S., MAGIC, and VERITA S. A firm detection 
of GRBs has not been reported so far llAlbert et al.l l2007l : 
lAharonian et al.1 l2009l : lAleksic et all |2010| ). but CTA may 
eventually achieve the goal with much lower energy thresh- 
old and much better sensitivity, which will also allow us to 
obtain time-resolved spectra of much better quality than 
ever before in the energy range greater than a few tens of 
GeV. In this work, we report the prospects for detecting 
long GRBs (durations longer than 2 sec) with CTA, and 
how the results depend on some of the array performance 
and GRB properties. The situation considered here is that 
GRBs are detected by some satellites and followed up by 
the LSTs, which will play a principal role in GRB obser- 



vations. In our simulation, we mainly assume Fermi/GBM 
as the burst trigger, which has already provided data with 
high statistics and may be active in the CTA era. We also 
make a rough estimate for the alerts from the SVOM satel- 
lite, which is planned to be launched before CTA operation. 
S ome complementary aspects of this work are also presented 
in llnoue et al] (12012 ). 

The organization of the paper is as follows. In sec- 
tion 2, we describe the intrinsic GRB properties of both 
prompt and afterglow emissions in our simulations. In sec- 
tion 3, we model the performance of the LSTs and the 
follow-up observation of GRBs. Performing Monte Carlo 
simulation, we estimate the detection rate of GRBs with 
the LSTs in section 4. Finally, section 5 is devoted to 
summary and discussion. We use cosmological parameters 
Ho = 70 km s"^Mpc"\ r2„ = 0.3 and = 0.7. 



2 INTRINSIC GRB PROPERTIES 

In this section, we describe our method of generating the 
intrinsic GRB samples in our Monte Carlo simulations and 
the assumed EBL model. It provides each sample with the 
prompt and afterglow properties using a luminosity function 
and well-known spectral correlations. We consider the GRB 
samples only with the duration > 2 sec and redshift < 5. 
It is also shown that our sample is consistent with some 
observational results of GBM and LAT onboard Fermi. 



2.1 Prompt emission 

We assume the luminosity function of GRBs, defined as the 
rate of GRBs per unit comoving volume at a redshift 2 and 
per logarithmic interval of peak luminosity Lp for I-IO'* keV 
in the cosmological rest frame, to be 

*(Lp,2)=p(2)0(Lp) . (1) 

Following IWanderman fc PiranI l|2010l ). we describe the func- 
tional forms of the GRB formation rate p{z) and the local 
luminosity function (pi^p) as broken power-laws: 



p{z) (X 



and 



(l){Lp) (X 



(1 + ^) 
(1 + ^) 



Z < Zp 
Z > Zp 



Lp < -L* 
Lp > L. 



(2) 



(3) 



In our fiducial case, 



S.ltni, OL = -0.171^;,, Pl = -lA4+°i and L. 



]^q52.5±o.2 gj.g g-1 g^^g adopted (jWanderman fc Piranllioiol ). 
The normalization is determined so that we reproduce the 
observed trigger rate for GBM of 250 yr~^ assuming its trig- 
ger condition as the peak photon flux of 1.5 ph s^^cm"^ in 
the 8-10^ keV band, satisfled by ~ 90 % of actual GBM 
bursts. In the following, for given values of Lp and z, we 
generate parameters of a burst regarding the prompt and 
afterglow emission. 

In order to model the prompt GRB observation with 
CTA, the simulated T90 distribution must be a good ap- 
proximation of the observed one, where Tgo is the time 
needed to accumulate from 5 % to 95 % of observed photon 
counts. The reason is that the delay time Tdday, which is 
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the start time of a follow-up observation as measured from 
the burst trigger time, is required to be shorter than the 
burst duration for the detection of the prompt emission. 
In order to determine Tgo of our simulated GRB samples, 
at first we assume that the duration of the prompt emis- 
sion in the energy range where the LSTs are most sensitive 
(from ~ 10 GeV to - 100 GeV) is equal to Tgo in the GBM 
band, and that the prompt phase provides constant gamma- 
ray luminosity Lave, evaluated by averaging the highly vari- 
able light curve. To determine the value Lave from the 
peak luminosity Lp, we take the average of ( - EiBo/Tgn)/L p 
in the observed events jGhirlanda et al.l l2004l . l2009l . l2010h . 
where Tgg = Tgo + z) and Siao is the isotropic-equivalent 
gamma-ray energy. The data of Tgo are taken from the GCN 
circulars. Then, we found (i5iso/Tgo)/Lp = 0.3±0.2. Despite 
the large scatter, we set 



0.3Lb 



(4) 



for all the simulated bursts. Next, we examine a co rrelation 
between log Lp and log Eiso (jGhirlanda et al.]|2012l ) employ- 
ing a chi-square fit with the effective variance weighting, and 
obtain 



log £150,52 = 1.1 log Lp,52 -f 0.56 



(5) 



with a standard deviation of 0.49 dex. Using this relation 
and taking the standard deviation into account, LJiao is de- 
termined by Lp in our simulation. Then, we determine the 
duration as 



^90 — L^iso/ La 



(6) 



and Tgo = il + z)T^o- 

As the intrinsic spectral shape of the GRB prompt 
e mission, we assume the Band function, which is given 
bv lBand et aP (|l993h as 



Ne oc 



L;'^exp(-£;/L;o) 



L; s; (a - l3)Eo 
E-^id- l3)Eo 



(7) 



where a and /3 are the high- and low-energy photon in- 
dices, respectively. If /3 < —2 and a > —2, then the pF^ 
spectrum has a peak at Ep — {2 + a)Eo. For given Lp, 
we determine -Ep accordin g to the Ep-Lp relation examined 
bvlGhirlanda et al.l (l2009h (the so called Yonetoku relation: 
lYonetoku et al.ll2004l ) taking the standard deviation into ac- 
count. The average luminosity Lave determines the normal- 
ization of the Band function. The photon indices a and /3 are 
determined according to the distrib ution that was actu ally 
observed for bright BATSE bursts (|Kaneko et al.ll2006t ). In 
our fiducial case, only the bursts with /3 < —2 are treated. 

In summary, the luminosity function provides Lp and z 
for each generated burst, which are subsequently related to 
Eiso, Lave, Tgo and Ep. Photon indices a and /3 are deter- 
mined independently. 

For a check of the validity of our method that generates 
GRB samples, we simulate the bursts triggering GBM and 
compare their Tgo and fluence distributions with observa- 
tions. The observed data (from GRB 080714 to 101130) were 
taken from the GCN circulars. In the top panel of Figure [U 
the scatter plot of the fluence and Tgo for both the simulated 
and the observed GBM samples are shown. In the bottom 
panel of the same figure, we show the simulated fluence dis- 
tribution for GBM bursts in the 8-10^ keV band superposed 
on the observed one. Both panels of Figure [T] show that the 
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Figure 1. Top panel: The scatter plot of fluence and Tgo for 
the simulated (black cross) and the observed (red square) GBM 
samples. Bottom panel: The simulated fluence S (8-10'^ keV) dis- 
tribution in the form of dAr/d(log S) for GBM bursts (black solid 
line) superposed on the observed one (red dotted line), which is 
normalized as Af = 1. 



actually observed properties of fluence and Tgo are well re- 
produced by our model. We find that 20 % of the GRBs in 
our simulation have Tgo larger than 100 sec, the typical value 
of Tdeiay assumed in our fiducial case as described later. This 
is only slightly larger than that for actually detected events 
(16 %). Our modelled duration distribution can be consid- 
ered to be in sufficiently good agreement with the observed 
one for our purposes of estimating the detection rates. 

Furthermore, we verify that our method roughly agrees 
with the Fermi/LAT detection rate, which is estimated 
as follows. First, we estimate the all-sky event rate to be 
~ 1000 yr~^ for the bursts satisfying Tgo > 2 sec and z < 5 
by deconvolving the observed GRB trigger rate of 250 yr~^ 
with the FOV of GBM (9.5 srj3 and the trigger efficiency 
of GBM obtained with our simulation. Next, we calculate 
the detection efficiency of LAT for the bursts satisfying both 
Tgo > 2 sec and z < 5 using our intrinsic GRB model, where 
a burst is judged to be detectable with LAT when the num- 
ber of detected photons > 100 MeV by LAT is larger than 
10. He re we take the LAT effective area from iAtwood et all 
(j2009l ) and consider the signal-dominated regime where the 
background can be neglected, and take into account its de- 
pendence on incidence angle simply by a factor of cos 6. 
The incidence angle is isotropically distributed in the LAT 
FOV. Then, combining the above derived values with the 



^ http:/ /gammaray.msfc. nasa.gov/gbm/instrumcnt/description 
/character. html. 
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Figure 2. Examples of the vFi, spectra of the prompt emission 
for the case of -Roxtra = (dotted curve) and 0.1 (soUd curve) are 
illustrated. For both cases, a = -1.0, /3 = -2.5, Ep = 250 keV, 
and 2 = 1.0 are assumed. Cutof fs near 100 GeV are c aused by 
the EBL attenuation modeled bv lRazzague et"aL I Hooi). 

LAT FOV of 2.4 sr (jAtwood et al.ll2009l ). we estimate the 
LAT detection rate to be 12 yr~^, which is roughly consis- 
tent with the observed one, about 7-8 yr~^ ((Granot et 'all 
I2OI0I ). Strictly speaking, our simulated LAT detection rate 
is slightly overestimated. Better agreement may possibly be 
achieved by accounting for spectral s oftening below the LAT 
band suggested for subset of events (|Beniamini et al.ll201ll ). 
more sophisticated model of LAT detection conditions, etc. 
We leave this issue as future work (though a brief discussion 
will be given in Section [S]). 

For the estimate of the detection rate with CTA/LSTs, 
we also simulate the case in which the prompt emission has 
an extra hard component in addition to the Band spec- 
tral component. Such an extra compo nent has been con- 
firmed by Fermi/LAT in some bursts (|Abdo et al.1 l2009bl : 
lAckermann et al]|2010l ). However, at present, the properties 
of the extra component, such as its spectral slope, amount 
of released energy, and the fraction of events with the ex- 
tra component, are highly uncertain. Therefore we introduce 
another parameter, iiextra = icxtra/iavc, where Lcxtra is the 
luminosity of the extra component in the 0.1-100 GeV band. 
Roughly speaking, iZoxtra is the ratio of the LAT band lu- 
minosity to t hat of GBM. For LAT bursts, iZextra is about 
a few 10 % (|Ghisellini et al.ll2010l ). and from the EGRET 
and LAT observations, the photon index in the GeV re- 
gion is about —2 on average (|Dinguslll995l :[ Ghisellini et al.l 
I2OI0I ). Below, for simplicity, we assume that all simulated 
bursts have the same 7?extra = 0.1 and a photon index of 
—2. Figure [2] illustrates examples of the vFi, spectrum of 
the prompt emission for the case of 7?extra ~ and 0.1. In 
both the cases, a — —1.0, 13 = —2.5, Ep = 250 keV, and 
z — 1.0 are adopted. Cutoffs n ear 100 GeV are caused by 
the EBL attenuation modeled bv lRazzaaue et al.l (|2009l ). We 
assume this EBL model in our fiducial case (see Section [27S)) . 

2.2 Afterglow emission 

High-energy afterglow emission in the CTA band is also 
highly uncertain. Then, foUowing lGhisellini et al.) (|2010l ). we 
assume a simple phenomenological model to describe the lu- 
minosity of the afterglow emission in the 0.1-100 GeV band 



in the cosmic rest frame as 

L.c(T')^10-ergs-(^)(|-)^' , (8) 

where pt is the temporal decay index, and T' — T/(l + z) 
is the elapsed time from the burst trigger in the cosmic rest 
frame. We choose the n ormalization to r eprod uce the re- 
sult seen in Figure 4 of iGhisellinT et al.l (|2010l '). Now the 
isotropic-equivalent energy and the burst duration of the 
prompt emission, Eiso and Tgo, are determined as described 
in the previous subsection. We assume Lag has nonzero 
value only if T > T90. In our fiducial case, we set the tempo- 
ral and spectral energy indices pt and pe (i.e., F^ oc T^ ' 1/''^ ) 
to —1.5 and —1.0, respectively l|Ghisellini et al.ll2010l ). 



2.3 EBL attenuation 

We need to take into account the attenuation of gamma- 
ray photons b y the EBL. In our fiduc ial case, we adopt the 
EBL model of lRazzaaue et al.l l|2009l ). This model is appli- 
cable for 2 < 5, so that only bursts with 2 < 5 are consid- 
ered in our simulation. As we see later in Figure (5] higher- 
redshift {z > 5) bursts are not expected to make a signifi- 
cant contribution to the detection rate with CTA (see also 
Ide Souza et al.ll201ll '). For comparison, we al so make calcu- 
lation s with the "Best Fit 06" EBL mode0 of lKneiske et al.1 
(|2004h at z < 5. In this model, the gamma-ray horizon, the 
location where the 77 optical depth r^.^, {E, z) = 1 lies at 
lower-redshift by a factor of ~ 1-2 at i5 < 100 GeV com- 
pared to our fiducial EBL model. However, we found that 
the resulting differences in the detection rate and the red- 
shift distribution for these two models is not large. 



3 MODELING OF GRB OBSERVATIONS 
WITH CTA/LSTS 

In this section, we model the CTA observation of GRBs to 
investigate the GRB detection rate and the expected GRB 
properties. The arrays for CTA will be constructed at two 
sites, each in the northern and southern hemispheres, and 
designed to have three types of telescopes: LSTs, MSTs, and 
SSTs. The LSTs will be the most crucial component for de- 
tecting G RBs because they domin ate sensitivity below 200- 
300 GeV (|CTA Consortium! [20T1I ') where EBL attenuation 
is expected to be less severe, and have the capability of fast 
slewing (180° in 20 sec). In this paper, we consider only 
one array site and the sensitivity of LSTs alone unless it is 
explicitly stated otherwise. 

We consider the situation where GRBs are detected by 
GBM and then are followed up by the LSTs, pointing to 
the centroid of the GBM error circle. Bursts that occur by 
chance in the LST FOV are not considered because their 
detection rate is expected to be very 



htt p : / / www .desy.de/" kneiske / downloads .html 

Assuming the all-sky GRB rate at the Earth of 10'^ yr~^, duty 
cycle of 10 %, the LST FOV of 4.6°, and the fraction of the 
detectable events in the FOV of 10 %, we get 0.004 yr"!. 
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3.1 GBM localization 

Follow-up observations are made only for sufficiently well- 
localized bursts. Therefore, we set a threshold value of the 
localization accuracy of the alerted burst for the follow-up 
observation. In the case of the MAGIC telescope, its FOV is 
3.5°, and t he threshold of the error ra dius to start the follow- 
up is 1.5° l|Garczarczvk et al]|2009l ). In our simulation for 
CTA/LSTs, we set the threshold value, oth, to 3.5° (i.e., 
the observation is made only when the alerted error radius 
with 1-sigma accuracy is less than (Jth)- The error radius of 
each sample is given as a function of energy fluence using 
the GBM burst data set, where the 1-sigma error radius of 
individual bursts is estimated as (a^tat +<'"aya)^''^ (= '^)- The 
systematic error ctbvb o f 3° is assumed for all the bursts for 
simplicity (|Briggs et al.. 2009). and the statistical error (Tstat 
is given by ground position data (V. Connaughton, private 
communication) which are reported > 1 min after the burst 
triggers. 

From these assumptions, the fraction of the triggered 
bursts for which the error radius is less than oth is — 
21 %. Note that in most cases, only the flight posi- 
tion data, which are more poorly determined compared to 
the ground one, are available during the prompt phase. 
Hence, as long as the current localization capability of 
GBM is considered, the simulated localization efficiency 
(i.e., 21 %) is larger than reality, which causes an overes- 
timate of the detection rate of the prompt emission. Nev- 
ertheless, we adopt this assumption expecting the improve- 
ment of the localization speed of GBM and its accuracy 
by the CTA era. For the afterglow phase, the above as- 
sumption for GBM localization has little influence. With- 
out any improvement of the GBM localization, nontrivial 
foUowup strategies that compensate for the limited FOV 
of the LSTs may be helpful to search for the prompt 
emission of G BM bursts, such as scanning the error cir- 
cle over time iFinnegan fc for the VERITAS Collaboration! 
I2OIII ). or divergent pointing of the LSTs over the error circle. 

The error radius of GBM localization is always larger 
than the radius corresponding to the LST's FOV of 2.3°. 
Hence in our simulation, we take into account the probabil- 
ity that the burst is in the FOV of LSTs after slewing for 
each burst alerted by the GBM. It is given by integrating 
the two-dimensional Gaussian distribution with the stan- 
dard deviation of cr from to 2.3°. In our fiducial case (i.e., 
o"th = 3.5°), the probability is about 20 %. The dependence 
of the detection rate on crth is shown in Section 14.31 



3.2 Delay time and detection conditions of CTA 

When a sufficiently well localized position of a GRB is 
obtained, the more rapidly CTA points toward the posi- 
tion, the more chances for GRB detection it has. In the 
case of MAGIC-I, if the distribution of the delay time be- 
tween the start of observations and the GRB trigger time 
(Tdoiay) is fitted by a log- normal probability density func- 
tion (PDF), then the PDF multiplied by Tdeiay has a peak 
at Tdeiay of ^ 160 scc and a standard deviation gd eiav of 
~ 0.5 dex jCarczarczvk et al.ll2009l : [Albert et"aIll2007^ . Tak- 
ing into account the ability of rapid slewing of the LST, we 
assume that Tdeiay of the LST obeys a log-normal distribu- 
tion with Tdeiay of 100 s and (Tdeiay of 0.4 dcx in our fiducial 
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Figure 3. The differential sensitivity curves of our CTA/LSTs 
model (red solid curve) and tlie official CTA (black dotted 
curve) jCTA Consortiumll2010l) . The official one includes the con- 
tributions from all types of telescopes (i.e., LSTs, MSTs, and 
SSTs), where at < 200-300 GeV LSTs dominate the sensitivity. 
The exposure time is 0.5 h. The right end of the horizontal axis 
(i.e., 300 GeV) is the maximum energy considered in our simula- 
tion. 



case. In addition we impose Tdeiay > 20 sec. It may be con- 
servative to use Tdeiay = 100 sec for LSTs considering that 
the average slewing time of the telescope (not the de lay time) 
of ~ 90 sec for MAGIC-I (|Garczarczvk et al.l [20091 ) is more 
than 70 sec longer than that for LSTs, < 20 sec. Simula- 
tion results for the other values of Tdeiay will be shown in 
Section 1321 

The last step for the GRB observations is to compare 
the gamma-ray flux with the sensitivity of CTA. Here we 
deflne N-^ as the total photon counts and A^'bg as the total 
background counts, where both are obtained by integration 
over the energy range from i?iow to 300 GeV for a given ob- 
servation time. The low-energy end of the integration i?iow 
is given in the next subsection. The upper bound of the in- 
tegration, 300 GeV, is sufficient for discussing the detection 
rate with LSTs. The exposure time is set to be Tgo — Tdoiay 
for the prompt emission and at most 4 hours for the after- 
glow. Fo llowing the simple est i mates of telescope sens itiv- 
ity (e.g.. lAharonian et al.ir200ll : ICTA Consortiumll201Gl ). we 
judge that a burst is detected if all of the following condi- 

2, and 



tions are satisfied: (1) A'^^ > A^niin, (2) N-y > 
(3) > eA^bg, where A'^min = 10, m = 5 and e = 0.05. 
The condition (1) concerns the minimum number of pho- 
tons required to create the sufficient Cherenkov radiation. 
The conditions (2) and (3) are about the statistical sig- 
nificance and the systematic error, respectively. To evalu- 
ate these conditions for each sample, we need the effective 
area of CTA/LSTs, background spectrum, angular resolu- 
tion, and -Eiow These are given in the next subsection. 



3.3 Performance of the LSTs 

In this subsection, we construct a toy model for the 
LSTs, which reproduces th e differential sensitivity of CTA 
that is publicly available (|CTA Consortium! l2010h . In or- 
der to calculate the differential sensitivity of our model, 
we set: the effective area of LSTs after all cuts for 
gamma rays A-y, the effective area for the background 
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Abg, the background spectrum Jbg in units of parti- 
cles s~^cm~^sr~^GeV~'^, and the angular resolution (f>. 
These are given as functions of photon or particle energy. 
The effective area is approximately the same as be- 
cause the background is dominated by electrons at energies 
below ~ 100 GeV, inducing electromagnetic showers simi- 
lar to gamma rays ||CTA Consortiumll2010l : lAharonian" et alj 

Then we assume ^bg oc A.y and Jbg c>c E 
according to the Fermi/LAT and the H ESS observations 
of cosmic-ray electrons and positrons (|Fan et al.l I2OI0I : 
lAharonian et alJ I2OO8I: lAbdo et alj|2009l). In addition, as- 
suming (j> oc 

g-0.5 (|CTA Consortiumll2010l ), we can cal- 
culate the differential count rate of the background as 
dRhg/dE = TrJbg^bgf?!'^ « E~'^-°^^Aj. The normalization 
of dRhg/dE is given below. 

Once the functional form of A-y is given, a shape of 
the differential sensitivity curve (that is, its energy depen- 
dence) is determined. We found that by using the effec- 
tive area of the MAGIC telesco pe at its trigger level as 
A^ (|MAGIC Collaboration! lioTlh H we can reproduce well 
the shape of the official sensitivity curve of CTA at < 200- 
300 GeV, where the LSTs dominate the sensitivity. There- 
fore as a functional form of at zenith angle Szen ~ 20°, 
we adopt the MAGIC effective area at the trigger level in 
our simulation. We normalize A^ introducing a factor of /a, 
where fA = l means that A^ corresponds to the effective 
area at the trigger level of MAGIC. 

For any value of /a, we can fit the sensitivity curve 
of our model to the public one with the normalization of 
dRhs/dE. In Figure [S] the official differential sensitivity of 
CTA is shown as the black dotted curve for the exposure 
time of 0.5 h, where contributions from all types of telescopes 
(i.e., LSTs, MSTs, and SSTs) are included, though LSTs 
dominate the sensitivity for E < 200-300 GeV. Because the 
sensitivity for this exposure time is limited by the statistical 
significance (see the detection condition (2) described in the 
last part of Section r3.2p at < 300 GeV, we can determine the 
normalization of dR^^^/dE as a function of /a. In the case 
of /a = 1, we need i?bg — 0.34 Hz for fitting, where _Rbg is 
calculated by the integration of dR^^g/dE from 20 GeV to 
20 TeV. The sensitivity curve of our model at O^en = 20° is 
shown as the red solid curve in Figure[3l One can see that the 
two curves are similar. This figure shows the photon energy 
up to 300 GeV since we calculate the photon counts below 
this energy in our simulation. Indeed, the GRB spectrum 
above ~ 100 GeV is expected to be severely attenuated by 
the EBL, so that our artificial cutoff at this energy does 
not affect the rate estimate. Above 200-300 GeV, the MSTs 
and the SSTs have better sensitivity compared to the LSTs, 
so that the sensitivity of our model deviates from the total 
array sensitivity. In order to keep the two sensitivity curves 
consistent with each other for arbitrary value of /a , we need 
Rbs = 0.34/a Hz to normalize dRhg/dE considering the 
detection condition (2). Although the value of /a influences 
the condition (1) and (3), it varies the expected detection 
rate only a little (see Section [4.4|) . In our fiducial case, we 
assume /a = 1- 

As 0zen increases, the shower maximum height gets 



^ The effective area at the trigger level decreases toward the low 
energy more gradually than after all cuts. 



higher and Cherenkov light is more absorbed, so that the 
area of light pool becomes large, while its density becomes 
low. In order to take into account the dependence of the 
effective area on the zenith angle, we use the following form: 



A-,{E,e. 



where E' 



E(cos 62., 



0' 



-A.y{E',0° 



(9) 



1.7 and C 



2.4 l|Aharonian et al]|l999l ). Also for the angular resolution 
0, we multiply by a factor of (-E'(e'zcn)/£'(20°))"^/^ 

The differential count rate of gamma rays becomes max- 
imum at ~ 60 (cos6lzon)"^ GeV in our model for a source 
with its photon index of —2, neglecting the EBL attenuation. 
We set -Eiow, which is the low-energy end of the integration 
for the photon counts, at an energy lower than this as 



-Elow {dz 



20 (cos^z 



GeV 



(10) 



As we describe in Section [4.41 the exponent of cosine, —3.3, 
has only a small infiuence on the detection rate. 



4 GRB DETECTION RATE WITH CTA 

In this section, we show the results of our Monte Carlo sim- 
ulation on the GRB detection rate with the LSTs, which are 
obtained under the assumptions described in Section [2] and 
SectionO Following is a summary of the simulation process. 
First, in order to generate intrinsic GRB samples, we ran- 
domly give Lp, z, a and /3 according to the distributions 
described in Section [2. II Using several correlations with re- 
spect to the prompt emission properties, we determine Ep, 
Eiso, iave, and Tgo from Lp; for the first two parameters, 
the deviation from the best fit line of the correlation is given 
by a Gaussian random variable. Then, we can calculate the 
prompt and afterglow gamma-ray fiux arriving at the Earth 
for each generated burst taking into account the EBL at- 
tenuation. Second, we set the trigger and the localization 
conditions to determine the detected and localized events. 
The former is given in terms of the peak photon fiux, while 
the latter is given by the probability of sufficient localization 
as a function of the fiuence (deduced in Section I3.1|l . Then 
it is judged whether each sample is localized sufficiently well 
and whether it is in the FOV of the LSTs after slewing to 
the best position. Finally, at the stage of follow-up observa- 
tion by LSTs, Tdeiay is given by a log-normal distribution, 
and Szcn is isotropically distributed independently of Tdday, 
where we assume a 10 % duty cycle and 6'zon < 60° as the 
observational criteria. We evaluate the detection conditions 
described in Section [3.21 on each sample taking into account 
the zenith angle dependence of the array performance (as 
shown in Section [3. 3p . 

In this simulation, we consider only one array site and 
the GBM alerts (with its trigger rate of 250 yr~^) alone. 
We limit our simulation to the GBM bursts with Tgo > 
2 sec and z < 5, the fraction of which is ~ 80 % of all 
GBM bursts. Note that this criterion does not greatly affect 
our final result of the detection rate. One of the reasons 
is that the trigger rate of short GRBs (T90 < 2 sec) and 
high-redshift GRBs {z > 5) is typically smaller than the 
total trigger rate. Second, it is nearly impossible to start 
foUowup observations of short GRBs within their duration 
of prompt emission since it takes more than 2 sec to receive 
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the positional information from satellites. In addition, short 
bursts have smaller fluences than the long bursts, which will 
make the localization of the short GRBs with GBM more 
difficult. Finally, detection rate of the high-redshift GRBs is 
expected to be smaller (especially for afterglows) due to the 
severe EBL attenuation of gamma rays. 

Hereafter, we classify generated GRB samples into sev- 
eral sets, in accordance with steps from the trigger by GBM 
to the detection by CTA: 

Alert : GRB samples that are detected by GBM and satisfy 
T90 > 2 sec and z < 5. 

CTAobs : Samples belonging to Alert whose positions are 
determined within zenith angle 6'zon < 60° with the error 
radius a smaller than (Tth. Moreover, the CTA duty cycle of 
10 % is also taken into account to select the samples of this 
class. 

Fobs : Samples belonging to CTAobs which satisfy the cri- 
terion Tdelay < Tbo • 

Pdet : Samples belonging to Fobs which are actually lo- 
cated within the FOV of LSTs and their prompt emissions 
are detectable by LSTs. 

Adet : Samples belonging to CTAobs which are actually 
located within the FOV of LSTs and their afterglows are 
detectable by LSTs. 

LSTs makes observations for the CTAobs bursts. 



4.1 Detection rate in the fiducial case 

First, we show results of the detection rate with LSTs for our 
fiducial case in which we assume: (1) only the Band com- 
ponent for the prompt spectrum, i.e., no extra hard com- 
ponent (7?oxtra = 0), (2) GRBs whose high-energy photon 
index of the Band component is /3 < —2, (3) for the after- 
glow emission, the temporal index pt of —1.5 and the spec- 
tral energy index o f —1, (4) the EBL model provided by 
iRazzague et al.l (|2009l '). (5) the typical delay time of start- 
ing observation, Tdeiay, of 100 sec and crdeiay of 0.4 dex, and 
(6) the threshold value of the error radius ath of 3.5° to 
start the follow-up observations. The dependence of the de- 
tection rate on parameters related to the LST performance 
and GRB emission is quantitatively discussed in later sub- 
sections. 



4-1.1 Frompt emission 

With the above assumptions, we obtain 0.03 yr~^ as the de- 
tection rate of the GRB prompt emission with LSTs (see the 
first line in Table [T]). Factors reducing the rate of Alert to 
that of Fdet are as follows. Among the samples belonging to 
Alert (~ 200 yr"^), 21 % (i.e., 43 yr"^) are localized enough 
to be followed- up by LSTs (see Section [3.11) . Moreover, im- 
posing ^zen < 60° and duty cycle of 0.1 we get 1.1 yr~^ for 
CTAobs sample. Of these. Fobs GRBs are 0.45 yr~^, and 
6 % of Fobs GRBs belong to Fdet. Note that only one ar- 
ray site and only the alerts from GBM are assumed in this 
calculation. We discuss more general cases in Section [5] 

In Figure m we show the cumulative distributions P{> 
N~f) of the total photon counts A'^., where A'^, is calculated by 
integration over -Eiow < E < 300 GeV. The red solid curve 
shows Fdet samples for the fiducial case. Once the prompt 
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Figure 4. The cumulative distributions of the total photon 
counts N^i, P{> Nj), where A''-, is calculated by integration over 
^'low < E < 300 GeV. The low-energy end of the integration Ei^^ 
is defined by Eq. I|10|l . The red solid curve shows Pdet samples for 
the fiducial case (see Section |4] or Table [T] for explanation of clas- 
sification for our samples). The blue dashed curve represents the 
case of -Roxtra = 0.1 for Pdet samples with the other parameters 
fixed to the fiducial values. 

emission is successfully observed, we can expect the photon 
counts A-y > 10^ with the probability of 60 %. We also calcu- 
lated the photon counts (> 1 GeV) expected for fermi/LAT, 
and compared with N-^. Then it is found that CTA/LSTs 
can detect 10-100 times more GeV photons than LAT for 
almost all the Fdet samples. Therefore CTA can provide the 
temporal and the spectral structure of high-energy emission 
of GRBs with higher significance than any other current in- 
struments. 

Figure [S] shows the redshift distribution of Pdet samples 
for the fiducial case with the red solid line in the form of the 
PDF. We found that for Pdet, z ~ 1.5 is typically expected 
and 90 % have redshifts less than 3.5. These results are little 
infiuenced by our limitation of z < 5, while it can depend on 
the assumed EBL model. Again we find that t he difference in 
the re dshift distrib ution between the mod el of lKneiske et al.l 
(2004') and that of'Raz zague et al.l l|2009 l) is very small (see 
also Table [2] with regard to the detection rate). For compar- 
ison, we show the result for CTAobs samples (that have not 
been selected by the CTA performance) with black dotted 
line. We can see that the distribution of CTAobs is close to 
that of Fdet, which implies that the redshift distribution of 
Fdet events is mainly determined by the GBM sensitivity. 

4.1.2 Afterglows 

The detection rate of the afterglow is obtained as 0.13 yr~^. 
Among CTAobs samples, 12 % satisfy detection criteria. De- 
tections of the afterglow are expected for ~ 84 % of GRBs 
whose prompt emission is detectable. Table [T] summarizes 
the results for our fiducial model assumptions. A similar es- 
timate for the SVOM satellite will be presented in Section 
1431 

In Figure [5] blue dash-dotted lines show the redshift 
distribution of the Adet samples. We can expect typically 
z ~ 1 for Adet. It is found that the peak of the distribution of 
Adet samples is shifted toward lower z from the distribution 
of CTAobs samples. This implies that in contrast to Pdet, 
low-z samples are selected by the CTA sensitivity. We found 
that 90 % of Adet samples have redshifts less than 2.9. 
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Table 1. Expected event rates for one array site. We separately consider the cases in which each of Fermi/GBM and SKOM/ECLAIRs 
are the alerting detectors. Our fiducial model parameters are assumed. Each class is defined as follows: {Alert) Samples of GRBs that are 
detected by GBM and satisfy T90 > 2 sec and 2 < 5; [CTAobs) GRB samples belonging to Alert whose positions are determined within 
zenith angle 9zon < 60° with an accuracy of cr < 3.5°. Moreover, we take into account CTA duty cycle of 10 % to select CTAobs samples; 
{Fobs) GRB samples belonging to CTAobs which satisfy the criterion T^^^i^y < T90; [Pdet) GRB samples whose prompt emissions are 
detectable by LSTs; {Adet) GRB samples whose afterglow emissions are detectable by LSTs. 

Alert CTAobs Fobs Fdet Adet 
Fermi/GBM [yr^i] ~ 200 LI 045 003 0.13 

SVOM/ECLAIRs [yr-^] ~ 56 2.0 0.65 0.1 0.37 
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Figure 5. The simulated redshift distributions are shown as the 
form of the probability density function (PDF). For Pdet samples 
(see Section|4]l, the red solid line describes the fiducial case, while 
the blue dash-dotted line is for the case of i?oxtra = 0.1 with the 
other parameters fixed to the fiducial values. The distribution of 
Adet samples and CTAobs samples in the fiducial case are rep- 
resented by the blue dash-dotted line and the black dotted line, 
respectively. 

4.2 Dependence on delay time distribution 

In Figure [G] we show the dependence of the GRB detection 
rate on the typical delay time, Tdciay, which is introduced in 
Section [3.21 The top and the bottom panels represent cases 
for the prompt emission and the afterglow, respectively. In 
both panels, the horizontal axis represents Tdeiay, whereas 
the vertical axis shows the ratio of the detection rate to 
that for the fiducial parameter set. Hence, in both panels, 
the curves labeled as fiducial have the ratio of 1 if Tdeiay = 
100 sec. 

First, let us consider the prompt emission (the top panel 
of Figure O. The red solid curve represents the fiducial case 
except for Tdeiay If Tdeiay = 60 sec for LSTs, the detection 
rate is enhanced by a factor of 1.3. The light-blue dashed 
curve shows the result for the case where the extra com- 
ponent with JJcxtra = 0.1 is added to a spectrum of the 
prompt emission in all GRB samples. In this case, indepen- 
dently of Tdeiay, the detection rate is doubled compared to 
our fiducial case (-Rextra = 0). In addition, to see the in- 
fiuence of the dispersion of Tdciay distribution we draw the 
magenta dotted curve for the extreme case of no dispersion, 
i.e., Tdeiay = Tdeiay for all eveuts. At Tdciay = 90 sec where 
it is comparable to the peak of the Tgo duration distribu- 
tion of CTAobs GRBs, the two lines cross each other. If 
Tdeiay 90 sec, the dispersion makes the events with Tdeiay 
smaller than the central value Tdeiay, which enhances the 




Tdeiay [sec] 

Figure 6. Dependence of the GRB detection rate on Tdeiay, 
where Tdday is the typical delay time between the follow-up ob- 
servation and the trigger. The right- and the left-hand side of the 
vertical axis represents the detection rate and its ratio to that for 
our fiducial case (i.e., Tdday = 100 sec), respectively. Top panel: 
Results for the prompt emission. The red solid curve represents 
the fiducial case except for Tdeiay ■ The magenta dotted curve rep- 
resents the case with no variance of the delay time distribution. 
The light-blue dashed curve represents the case in which the ex- 
tra spectral component with i?extra = 0.1 is introduced. Bottom 
panel: Results for the afterglow emission. The blue solid curve rep- 
resents the fiducial case except for Tdeiay The brown dot-dashed 
and the pink dashed curves represent the cases in which the af- 
terglow temporal index pt is set to —1.3 and —1.8, respectively. 



detection rate. On the other hand, if Tdciay lS 90 sec, the 
dispersion makes the events with Tieiay larger than the cen- 
tral value Tdciay, which reduces the detection rate. 

Next, let us consider the afterglow (the bottom panel of 
Figure |6} . As the red solid curve in the top panel, the blue 
solid curve is for the fiducial case except for Tdeiay This 
curve is for the afterglow temporal index pt of —1.5, while 
the brown dot-dashed and the pink dashed curves are for 
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pt = —1.3 and —1.8, respectively. These suggest that for 
the afterglow detection, the important factor is not Tdeiay 
but others such as the well-localized alert rate and the low- 
energy sensitivity. 

The detection rate of the prompt emission is more sen- 
sitive to Tdeiay than that of the afterglow as shown in Fig- 
ure [6l This simply comes from the fact that Tdeiay affects 
the number of Fobs bursts that satisfy Tdeiay < Tgo. To see 
this explicitly, we show in Figure [7] the Tdeiay distributions 
of CTAobs, Fobs, Fdet and Adet in our fiducial case. We can 
see that the Tdeiay distribution of Fdet is about the same as 
that of Fobs, and they are shifted toward the shorter Tdeiay 
from that of CTAobs. On the other hand, the distribution 
of Adet is only slightly off set from that of CTAobs. In the 
afterglow case the condition Tdeiay < Tgo is not required. 

Although shortening Tdeiay to less than 100 sec does not 
have much impact on the improvement of the detection rate, 
we should note that shorter delay time is necessary for de- 
tecting the prompt emission. Because the minimum cutoff 
of Tdeiay is fixed to be 20 sec in our simulation, even when 
larger Tdeiay is assumed, shorter Tdeiay plays a significant role 
as mentioned above in the explanation of (Tdeiay = 0. Actu- 
ally, in the bottom panel of Figure[71 for Tdeiay = 100 sec, we 
can see that Tdeiay is distributed at less than ~ 100 sec in 
most cases where the prompt emission is detectable. There- 
fore it is crucial to keep the delay time as short as possible if 
the duration of the bursts in the CTA band and in the GBM 
band are similar to each other. In order to reduce Tdeiay, fast 
alerts with good localization is equally important to a rapid 
slewing of LSTs. As an example, we simulate the case in 
which Tdeiay cau not be shorter than 100 sec with the other 
parameters fixed to the fiducial values, and find that the de- 
tection rate of the prompt emission decreases by a factor of 
2 in this case. 
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Figure 7. The delay time distributions for each class in the 
case of our fiducial parameter set (see Section |4] or Table [T] for 
explanation of classification for our samples). The bottom panel 
and the top panel are in the form of the cumulative distributions 
and diV/d(log 

Tdeiay) normalized as — 1, respectively. Each 
line represents the different class: for CTAobs, the black dotted 
line; for Fobs, the black dashed line; for Pdet, the red solid line; 
and for Adet, the blue dot-dashed line. 



4.3 Dependence on the criterion of the error 
radius to start follow-up observations 

One of the criteria for LSTs to start follow-up observations 
concerns ath, whose fiducial value is set to 3.5° (see Sec- 
tion 13. 1[) . Most of the GBM alerts have the error radius 
larger than this fiducial value. The larger ath we use, the 
more chances of the follow-ups we have, while the efficiency 
of the detection decreases. This is because the bursts with 
larger error radii are less probable to lie in the FOV of the 
LSTs. 

In Figure [S] we show the dependence of the detection 
rate on ath. The results are shown as the ratio to the fidu- 
cial case. The red solid line (Fdet) and the blue dot-dashed 
line (Adet) show the increase in the detection rate by a fac- 
tor of 1.2-1.3 for (Jth = 5° compared to 3.5°. On the other 
hand, one should keep in mind that the more rapid increase 
of CTAobs (black dotted line) shows the decline in the de- 
tection efficiency (i.e., the ratio of Fdet or Adet to CTAobs) 
by a factor of 0.6-0.7. On the contrary, if we take ath to be 
less than 3.2°, most of the GBM alerts do not satisfy this 
criterion, so that the detection rate decreases very rapidly. 

It can be seen from Figure [8] that the detection rates of 
Fdet and Adet saturate around ath = 5°, so that it seems 
better to set ath to near this value as long as we take a 
simple strategy of the follow-up observations, i.e. all 4 LSTs 
point toward the centroid of the GBM error circle. 




3 4 5 

Threshold of error radius (Jth [deg] 

Figure 8. Dependence of the event rate on the threshold error 
radius ath, which is one of the criteria for LSTs to start the follow- 
up observations introduced in Section lij.il The results are shown 
as the ratio to our fiducial case. The red solid line represents the 
result for Pdet, the blue dot-dashed line for Adet, and the black 
dotted line for CTAobs (see Section |4] or Table [T] for explanation 
of classification for our samples). 
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Table 2. Parameter dependence of the detection rate for one array site of CTA/LSTs, where only GBM is assumed as the alerting 
detector. The cases other than the fiducial one are for those in which one of the model parameters takes different value from the fiducial 
parameter set (with the other parameters fixed). Note that the localization accuracy of GBM is estimated optimistically for the prompt 
emission phase compared to the current localization accuracy. See Section 14.41 for the explanations of each Case. 

fiducial Case 1 Case 2 Case 3 Case 4 Case 5 Case 6 Case 7 Case 8 

Tdclay = 60 s /3 > -2 -Rcxtra = 0.1 Knciske pf = -1.3 pt = -1.8 = -0.5 = -1.5 

Pdet [yr-i] 003 004 O04 O06 O02 

Adet [yr-i] 0.13 0.14 - - 0.11 0.15 0.10 0.19 0.03 
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. ^ A ^ 







Case 



Figure 9. The results in Table [2] are plotted as the ratio of the 
detection rate to the fiducial case. Each case corresponds to those 
of the table [2] The red solid and blue dot-dashed lines represent 
the results for Pdet and Adet, respectively. The same plot but for 
the simulated LAT detection rate, which is expected as 12 yr~^ 
in our fiducial case, are shown with the filled circle for Case 2, 3, 
and 4. 



4.4 Dependence on other parameters 

Table [2] summarizes how the detection rate changes when 
one of our model parameters takes different values from that 
in our fiducial parameter set (with the other parameters 
fixed). Each column describes the result for different cases: 
{first column) The fiducial case. (Case 1) The typical delay 
time of the follow-up rdday is set to 60 sec. (Case 2) The 
case including /3 > —2 according to the BATSE observa- 
tion. (Case 3) The extra spectral component of Rcxtrm = 0.1 
is assumed for all th e samples. (Case 4) The EBL model by 
iKneiske et al.l (|2004l ') is used. (Case 5 and Case 6) The tem- 
poral index of the afterglow pt is taken as —1.3 and —1.8, 
respectively. (Case 7 and Case 8) The spectral energy index 
of the afterglow pe is taken as —0.5 and —1.5, respectively. 
Figure [9] shows the same results presented in Table [2] but 
plotted as the ratio to the fiducial case. The red solid and 
blue dot-dashed lines represent the result for Pdet and Adet, 
respectively. 

The detection rate for Case 3 is about twice as large as 
the fiducial case. This fact remains unchanged for different 
Tdclay, which is represented in the top panel of Figure[6]with 
the light-blue dashed curve. This is because -Rextra of 0.1 
makes the bursts with intrinsically soft Band spectrum with 
/3 < —2.5 detectable; these soft bursts account for about 
half of Alert and CTAobs samples. By contrast, i^oxtra = 0.1 
has little infiuence on the photon-count distribution and the 
redshift distribution. Comparing Case 2 with Case 3, the 



latter shows a larger increase of the detection rate compared 
to the former. This simply reflects that the fraction of bursts 
with > —2 is less than those with /3 < —2.5. 

We also calculate the expected LAT detection rate of 
Case 2 and Case 3, and then obtain 24 yr~^ and 14 yr~^, 
respectively. The former is more than three times as large as 
the observed LAT rate for the bursts of Tgo > 2 sec and z < 5 
(about 7-8 yr~^), which seems unrealistic. The calculated 
LAT event rate for our fiducial case and Case 3 are similar to 
each other and about 1.5-2.0 times as large as the observed 
one. These somewha t large rates appear to be consistent 
with the analysis by iBeniamini et al.l (|201ll ). claiming the 
existence of the spectral softening below the LAT band in 
some of bright bursts, since we have not taken into account 
any spectral softening feature in our simulations. We discuss 
this point a bit more in Section [S] 

In addition to the cases summarized in Table[2] we stud- 
ied the dependence of the detection rate on the luminosity 
function (see Eq. ^ and (O), the normalization of the ef- 
fective area, /a (see Section |3.3|) . and the low-energy end 
of the LSTs sensitivity, Siow (see Eg. llOf) . For the luminos- 
ity function, we found that when one of the 6 parameters 
included is varied from the best fit value within its errors, 
the rate changes at most by a factor of ~ 0.8-1.3 for both 
Pdet and Adet. For /a, we see its infiuence on the detec- 
tion rate by varying the value from 0.3 to 3. At this time, 
as a function of /a, the background count rate 7?bg shifts 
from about 0.03 Hz to 3 Hz to keep the sensitivity of our 
LSTs model consistent with the official one. We found that 
this range of /a changes the detection rate by a factor of 
0.8-1 for Pdet, though the change of the Adet rate is neg- 
ligible. This refiects the fact that for non-detected prompt 
emissions, the detection conditions (1) and (2) (which are 
described in the last part of Section [3.2[) . have comparable 
importance, while for non-detected afterglows, the condition 
(2) is the most strict. Note that the expected photon counts 
7V^ can vary in proportion to /a. For iSiow, we see the detec- 
tion rate for the various exponent of cosine in Ea. (|10p . By 
varying it from —4.3 to 0, we found that the detection rate 
changes by a factor of ~ 1.0-1.2 for Pdet and ~ 0.9-1.1 for 
Adet. Therefore the luminosity function, /a, or Siow is not 
more sensitive to the detection rate than the other source 
parameters. 



4.5 CTA detection rate for SVOM GRBs 

Finally, we roughly estimate the detection rate for the case 
of alerts from SKOM/ECLAIRs, whose launch befor e the 
CTA operation is being planned (|Schanne et al.|[2010l ) and 
which wi ll provide well-lo calized alerts (< 10') of about 
80 yr~^ (|Paul et al.ll201ll ). Here we assume that its dura- 
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tion distribution is the same as the Burst Alert Telescope 
(BAT) onboard Swift because the t wo detectors, ECLAIR s 
(4-250 keV) and BAT (15-150 keV: ISakamoto et alJlioOTl ). 
cover a similar energy range. The Swift data on Tgo 
was taken from the web site (from GRB 041217 to 
GRB 110519A). 80 % of ECLAIRs GRBs are expected at 
z < 6 l|Paul et al.ll201ll ) and ~ 90 % of the BAT bursts are 
long GRBs. Hence we assume that the fraction of all GRBs 
triggering ECLAIRS that have Tgo > 2 sec and 2; < 5 is 
70 %. We set the delay time Tdday of 80 sec for all bursts 
for simplicity, while we assumed typically Tdeiay = 100 sec 
for GBM. This is because the SVOM alerts are expected to 
befaster than the GBM ones (i.e., < 1 min (jSchanne et al.l 

[ioS)). 

With the above assumptions, we can estimate the frac- 
tion of well localized events for which Tgo > Tdday is ~ 33 % 
for long GRBs 0. It is qualitatively expected that the frac- 
tion of Pdet or Adet to the bursts in the LST FOV localized 
by ECLAIRs is lower than that by GBM. The reasons for 
this are as follows: (1) ECLAIRs does not need high fluence 
for localization as much as GBM, which leads to a lot of 
burst samples with dim flux in the CTA band. (2) ECLAIRs 
has better sensitivity for softer bursts rather than the hard 
ones, which again leads to the small number of events de- 
tectable with CTA. Also, the redshift distribution of bursts 
localized by ECLAIRs shifts to higher redshift than that 
by GBM, so that larger fraction of events are severely af- 
fected by the EBL attenuation. In evaluating the fraction of 
Pdet or Adet to the bursts in LST FOV quantitatively, we do 
not simulate in detail the follow-up observation of ECLAIRs 
bursts with CTA taking into account its performance, but 
we just assume that the efficiency is the same as the case of 
Swift/3AT - if BAT is assumed to be the burst trigger and 
its trigger threshold for peak photon flux is simply set to 
0.4 ph s"^cm"^ in the 15-150 keV band (about 90 % of ac- 
tually detected BAT bursts have their peak flux above this 
value), we obtain lower detection efficiency than the case of 
GBM by a factor of 0.54 for the prompt emission and 0.34 
for the afterglow in the case of ath = 3.5°, Tdeiay = 80 sec 
for BAT, and 100 sec for GBM. Additionally, we must con- 
sider that SVOM operations feature a bias of preferentially 
pointing toward the anti-solar direction, which increases 
the probability of follow-up observations with ground-based 
telescopes. We set the i ncrease of the detect ion rate by a 
factor of 1.4 referring to lCilmore et al.] l|2010l ). 

With the above assumptions we can estimate the CTA 
detection rate for alerts from ECLAIRs as ~ 0.1 yr~^ for 
the prompt emission, while ~ 0.37 yr^^ for the afterglow 
(see Table [T]), which is about three times higher rate than 
GBM. Hence SVOM will probably become better for the 
detection of GRBs with CTA than Fermi. Note that the 
GBM localization accuracy is estimated optimistically for 
the prompt emission phase (see Section [3.ip . and also note 
that, unless the SVOM alerts are transmitted to enough 



5 http: / /swift. gsfc.nasa.gov/docs/swift/swiftsc. html 
^ GBM requires higher fluence to localize bursts, so that the du- 
ration of the localized bursts tends to be longer than that of the 
detected bursts to some extent. This is responsible for smaller 
ratio of bursts with Tgo > Tdeiay to the localized bursts for BAT 
(i.e., 33 %) than that for GBM (i.e., 43 %). 



ground stations, the delay time of 80 sec for SVOM alerts 
turns out to be an optimistic assumption. 



5 SUMMARY AND DISCUSSIONS 

In this paper, we have presented the prospects for detec- 
tion of GRBs with the LSTs of CTA, the most vital com- 
ponent for GRB observations with their fast slewing capa- 
bility and the best sensitivity at the lowest energies. We 
have modeled and simulated the follow-up observation of 
GRBs with T90 > 2 sec and z < 5 alerted by Fermt/GBM. 
Our GRB population is modelled according to a given lu- 
minosity function that is consistent with Swift observations 
together with well-known spectral correlations. We note the 
following strengths of our model: (1) It reproduces the flu- 
ence and duration distributions of GBM bursts (Figure [T]). 
(2) The differential sensitivity of our m odel LSTs is con- 
sisten t with the official one given by the ICTA Consortium! 
(|2010l ) in the energy range less than a few 100 GeV (Fig- 
ure [3]). (3) The Fermi/L AT detection rate is predicted to be 
12 yr~^ for the bursts satisfying both Tgo > 2 sec and z < 5, 
which is roughly consistent with the actual rate of 7-8 yr~^. 

Assuming Fermi/GUM alerts alone, our fiducial param- 
eter set predicts the GRB detection rate with LSTs for one 
array site (i.e., north or south) as 0.03 yr~^ for the prompt 
emission and 0.13 yr~^ for the afterglow emission (Table [T]). 
The expected event rates become larger when two array 
sites of CTA and additional alerts from SVOM/ECLAIRs 
are taken into account — if the array performance for the 
two sites are the same in the energy range less than a few 
100 GeV, the detection rates go up to about 0.3 yr~^ and 
1 yr~^ for the prompt and afterglow emissions, respectively, 
where we assume no overlap of FOVs of GBM and ECLAIRs. 
Note that in the above estimates, we optimistically esti- 
mate the onboard-localization ability of GBM in the prompt 
phase, expecting improvement by the CTA era (see Section 
13. 1|) . For the afterglow emission, our treatment of the GBM 
localization ability hardly affects the rate estimation. 

For our fiducial assumptions, once CTA succeeds in de- 
tecting the prompt emissions for GBM alerts, the total pho- 
ton counts A^.^ are expected to be > 10^ for 60 % of CTA 
detected events. Our simulation also shows that A'^, is 10- 
100 times larger than the number of GeV photons expected 
by LAT for the same bursts. This suggests that CTA can ob- 
tain the time resolved light curve in the energy range greater 
than a few tens of GeV with higher statistics than ever be- 
fore. We expect 90 % of the prompt burst detected by CTA 
to have redshifts less than 3.5, and 90 % of the afterglows to 
have less than 2.9. Because of the follow-up observation after 
the GBM alerts, the redshift distribution is more affected by 
the GBM sensitivity rather than that of CTA. Hence, more 
frequent detections of high-z prompt emission are expected 
for SVOM alerts. 

Studying the dependence of the detection rate on our 
model parameters, we found the following results: 

(i) For the prompt emission, if all GRB samples have the 
extra power-law component with luminosity 10 % of the 
Band component, the detection rate with CTA increases by 
a factor of 1.9. 

(ii) For the afterglow, the spectral index has a relatively 
large effect on the detection rate. It decreases by a factor 
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of 3 with the extreme assumption that all bursts have the 
softer spectral energy index of —1.5, while in the harder case 
of —0.5, the rate increases by a factor of 1.5. 

(iii) If the typical value of the delay time between the start 
of observations and the GRB trigger time is 1 min, the de- 
tection rate of the prompt emission increases by a factor of 
about 1.3 compared to our fiducial case. On the other hand, 
the rate of the afterglow depends slightly on the delay time. 
Although the shorter delay time does not have much effect 
on the enhancement of the detection rate itself, the ability of 
fast slewing (< 100 sec) is necessary for catching the prompt 
emission. 

(iv) It seems better to make follow-up observations for the 
alerts with the error radius of up to ~ 5°, as long as all 
4 LSTs simply point toward the best position localized by 
GBM, although the detection efficiency is decreased. 

(v) The detection rate varies by less than 30 % for variation 
of the other parameters including those of the luminosity 
function. 

Note that the LAT detection rate is estimated to be 
14 yr~^ in the case of (i), which is about 1.8-2.0 times 
as large as the observed LAT event rate. This discrepancy 
may be due to the existence of the spectral softening be- 
low the LAT band, s uggested for subset of bright events 
l|Beniamini et al.]|201ll ) , our crude (but reasonable) LAT de- 
tection conditions, etc. Even if the spectral softening effect 
is a dominant cause, the detection rate with CTA is not 
necessarily expected to decrease from our results by a simi- 
lar factor, 1-2, due to the following possible reasons. First, 
the spectral softening below the LA T band have been onl y 
suggested for some of bright bursts (jBeniamini et al]|201ll l. 
CTA will be sensitive even for typical (less bright) bursts, 
which might not have the spectral softening. Second, the 
bursts with the steep Band component below the LAT band 
may be hardly detected by LAT, whereas they might have an 
extra hard component at > 10 GeV, which can be relevant 
for CTA. We have not taken into account such possibilities 
in our simulations, which remain as future work. In any case, 
CTA may clarify the properties of the extra hard component 
and spectral softening of the Band component, which would 
strongly constrain the prompt emission mechanism. 

It is crucial to increase the alert frequency with good 
localization in order to achieve higher detection rates. Hence 
the improvement of the localization accuracy by future GRB 
alert facilities is important. Let us consider the alerting de- 
tector which has the same characteristics (e.g., sensitivity) 
as GBM except for the localization ability, and simply as- 
sume that the alerts from the detector have a constant er- 
ror radius a for all the alerted bursts. In this case, if we 
can use the localization error o less than the threshold 
ath, the probability of localization in the FOV gets higher, 
while the frequency of follow-up observations do not vary 
(CrAo6s~ 5yr-^).The increase of the probability saturates 
at cr ~ 1°, and then the detection rate of Pdet and Adet for 
this detector are enhanced to ~ 0.2 yr~^ and ~ 1.2 yr~^, 
respectively. It implies that future alerting detectors with 
the better localization ability are more desirable. 

As countermeasures against the alerts with large local- 
ization error, two observing strategies are proposed: one is 
to point each LST at different directions in the GBM error 
circle, while the other is to scan the region instead of point- 



ing one location. The optimal strategy for each observation 
can be determined by the estimate of the detection proba- 
bility, which requires the information on the trigger (such as 
the degree of localization error, the brightness, the expected 
delay time, and the zenith angle). In order to perform this 
and to increase the chance of detection, we need more de- 
tailed study on the sensitivity for each strategy as a function 
of some important parameters (such as zenith angle, the ex- 
tended FOV, scan speed, etc.), which remains as a future 
work. 

In the VHE region, the future Hig h Altitude Water 
Cherenkov (HAWC) Observatory mission (|Abevsekara et al] 
I2OIII I may also detect GRBs with its large FOV (~ 15 % 
of all the sky) and high duty cycle (~ 100 %), while the 
GRB spectrum may be difficult to determine by HAWC, be- 
cause of its small effective area and low energy resolution at 
~ 30 GeV compared to CTA (Goodman ct al. 2011). Hence, 
CTA and HAWC are complementary and both types of ob- 
servations are important^ 

Very r ec ently, iGilmore et all ()2011al lbh and 
iBouvier et al.l l|201lh independently studied on expec- 
tations for GRB detection rates by CTA, and we find their 
results are broadly in agreement with our conclusions here. 
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